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Abstract. Phonons, the quantum mechanical representation of lattice vibrations, 
and their coupling to the electronic degrees of freedom are important for understanding 
thermal and electric properties of materials. For the first time, phonons have been 
measured using resonant inelastic x-ray scattering (RIXS) across the Cu if -edge in 
cupric oxide (CuO). Analyzing these spectra using an ultra-short core-hole lifetime 
approximation and exact diagonalization techniques, we can explain the essential 
inelastic features. The relative spectral intensities are related to the electron-phonon 
coupling strengths. 
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1. Introduction 

It is generally accepted that the retardation of the Coulomb interaction due to the 
coupling between electrons and phonons is the mechanism for conventional BCS 
superconductivity. On the other hand, despite over two decades of unceasing research 
efforts, the mechanism for high-T c superconductivity in cuprates remains evasive. 
Although it has been suggested that electron-phonon coupling could be a key towards 
understanding this intriguing phenomenon (TJ [2} [3j 0], the topic is still being heavily 
debated [3 |8]. Whereas information about the phonon states can be obtained 
through inelastic neutron and x-ray scattering or Raman spectroscopy, the study of the 
coupling between electrons and phonons is more elusive. A number of experimental 
techniques provide an indirect measure of the electron-phonon coupling. For example, 
in angle-resolved photo-emission spectra of high-T c cuprates, kinks observed in the 
dispersion of the electron bands have been attributed to strong electron-phonon 
coupling. However, since the observed spectral features are dominated by electronic 
excitations, it is difficult to unambiguously assign the origin of these kinks to phonons 
rather than, for example, spin excitations [HI El [7]. Point-contact spectroscopy [10J, 
on the other hand, is sensitive to interaction of electrons with elementary excitations 
including phonons; however, this technique is not momentum resolved. In addition, 
the spectra strongly depend on the quality of the tunnel barrier, transmission matrix 
elements, and inelastic scattering in the barrier region. Resonant Inelastic X-ray 
Scattering (RIXS) has been shown to be sensitive to phonon excitations, but until 
now they have only been observed in combination with electronic excitations [9]. 

Here we demonstrate that RIXS can probe pure phonon excitations through 
the core-electron excitation, that provide element-specific and momentum-resolved 
information on the coupling between phonons and electrons, but that separates the 
phonon and the electron degrees of freedom in the final states. To achieve this, we 
employed RIXS at the transition-metal i^-edge. The underlying physical mechanism 
is as follows. The incoming x-ray excites an electron from the deep-lying Is core 
level into the valence band consisting of 3d and 4p states through quadrupolar and 
dipolar transitions, respectively. At the site where the absorption takes place, this 
excitation creates a sudden change in charge density, to which the lattice responds 
through the electron-phonon coupling. Since the deexcitation removes the electronic 
excitation created in the absorption, final states with only phonon excitations can be 
reached. Additionally, the very short core-hole lifetime prevents a complete relaxation 
of the lattice within the time scale of the resonant scattering process. 

2. Experiment 

In order to demonstrate the technique, cupric oxide (CuO) was used as a model system. 
CuO is of particular importance since it is the simplest member of the family that shares 
the same integral CUO4 plackets with the high-T c superconducting cuprates [TT]. The 



Observation of phonons with resonant inelastic x-ray scattering 



3 



CuO single crystal was grown and oriented along the (100) direction [T2] and measured 
at the 3-ID beamline of the Advanced Photon Source using a RIXS spectrometer based 
on a sapphire crystal analyzer at a back-scattering geometry [13]. The monochromatic 
beam of 22 meV resolution was obtained by the four-bounce asymmetrically cut Si (4 
4 4) monochromator. The sample was in vacuum to reduce the background due to air 
scattering, and the measurements were performed at room temperature. The overall 
energy resolution of the spectrometer was 38 meV (figure [1]) around the Cu K-edge as 
the test measurements of this new spectrometer reported in Ref. [IB]. The resolution 
measurements were repeated for each incident energy before and after the inelastic 
measurements to assure the observed effect was not due to a possible instrumental 
glitch. Spectra were collected at four incident photon energies ranging from the Is — > 3d 
(quadruple transition) to Is — > 4p (dipole transition). The incident energy extrema of 
8981 eV and 8997 eV were the limits of the instrument as the analyzer was designed 
to operate at a fixed Bragg angle of around 89.8° (very close to back-scattering), and 
analyzer energy was changed by thermal expansion of the scattering crystal, not the 
Bragg angle. The temperature of the analyzer crystal, ranging from 100 K to 500 K 
defines the energy, ranging from 8997 eV to 8981 eV, respectively. 

3. Results and Discussion 

Since there is evidence that electron-lattice coupling is strongest at the zone boundary 
in cuprates [14J, the measurements were taken in the first Brillouin zone such that the 
momentum transfer vector was along the (100) direction. This geometry corresponds 
to a scattering angle 26 of 16.93° with momentum transfer of 13.4 nm _1 (~1.8 nm _1 q 
resolution). Additionally, the phonon branches under investigation have their maximum 
energies at the zone boundary as reported earlier [15j . which makes them relatively easier 
to detect with the current instrumental resolution. 

The RIXS data are shown in figure [2] for several incoming photon energies. In the 
pre-edge region (8981 eV), which is dominated by Is — > 3d excitations, we observe a 
clear asymmetry in the spectral features with a maximum intensity around an energy 
loss of 45 meV. For larger incoming x-ray energies, the asymmetry remains but the 
intensity decreases. At an incoming x-ray energy of 8997 eV, the maximum coincides 
with the zero- loss peak with some asymmetry still present. It is known that the elastic 
signal also enhances as a result of ii'-edge resonance and one can argue that the observed 
spectral change may be due to this resonant behaviour of the zero-loss peak. However, 
it was not possible to fit the spectra according to this explanation. The effect of photon 
absorption coefficients cannot explain the reported anomaly either. The energy scan 
range is ±200 meV and the inelastic signal is observed within ±100 meV. The absorption 
coefficients are virtually the same for the short energy scanning range compared to the 
width of the absorption edge of ~20 eV. Therefore, the change in the spectra cannot 
solely be associated with energy dependence of the absorption coefficients, since photon 
absorption would only affect the overall intensity, not the line shape. 
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It is evident that the inelastic features show a clear resonant behaviour. In RIXS 
measurements, it is a common practice to take a spectrum away from the transition 
metal edge to prove the resonant behaviour [161 E] • However, in this case where the 
non-resonant measurements away from the Cu .fT-edge are not possible (see above for 
the discussion about the instrument and ref. p2]), the noticeable change in the spectra 
as a function of incident photon energy should be accepted as the necessary proof. 

There are several possible explanations for the origin of these inelastic features. We 
can readily discount crystal-field excitations, since the lowest transition is expected to 
be on the order of 1-1.5 eV [18]. Additionally, K-edge RIXS is known to be sensitive to 
magnon excitations [191 ED]- Since the RIXS process conserves spin, only two-magnon 
excitations are allowed with an expected combined energy of 0.3-0.5 eV [19l|2Q], which 
is an order of magnitude larger than the features in the present RIXS data. Therefore, 
we attribute these loss features to the excitations of lattice vibrations. The closeness 
of the energy loss of the RIXS spectra to the experimentally observed phonon modes 
supports our assertion [T5l [21] . 

As a result, it can be safely deduced that the observed variation in the spectral shape 
as a function of incident photon energy can be explained by the change in the relative 
intensities of the phonon modes against each other. At this geometry, three phonon 
modes are detectable with energies of 24 meV, 41 meV, and 70 meV [15]. We repeated 
these measurements along the (100) direction to make sure which branches are allowed 
for this geometry. The rough measurements were done using a 2.2 meV-resolution 
spectrometer that operates at 21 keV [22]. The experimental data demonstrate the 
coupling between the lattice and the electronic excitations created by the absorption 
process and this coupling leads to modulation in the phonon intensities. The electron- 
phonon interaction can be described by 



where oJ^ m creates a phonon with mode m, wavevector k, and energy fiu^m] after the 
absorption process, the coupling strength to the transient change in charge density 
with respect to the ground state is determined by Ak m . In the ground and final states 
this coupling is zero. This Hamiltonian corresponds to a displaced oscillator, and its 
solutions are well known. The sudden change in charge density due to the absorption 
process causes the lattice to respond through the excitation of multiple phonons. To 
understand the underlying physics, a single dispersionless optical phonon mode with 
energy hcoo and coupling A can be taken as an example. In the absorption process, 
phonons of this mode are excited following a Poisson distribution with the maximum 
given by A/hu , see the inset in Fig. [3j For a typical coupling A on the order of 100 
me V, A / Huq > 1 and it is expected that several phonons are excited in the intermediate 
state. One therefore expects that several phonons should be visible in the final state. To 
verify this, the RIXS cross section is calculated using the Kramers-Heisenberg equation 
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where \n), and |/) are the initial, intermediate, and final states, respectively; huo 
and Hoj' are the energies of the incoming and outgoing x-rays; T ~ p • A is the transition 
operator; and F is the full width at half maximum of the lifetime broadening of the core- 
hole. Figure [3] shows that through the Kramers- Heisenberg equation multiple phonon 
excitations are expected (see the calculations for F/huo = 3 and 5). Contrary to 
the calculations, the experimental spectra do not show multiple phonon excitations. 
This peculiar effect can be explained through numerical calculations by increasing 
T/huo. Even though the intermediate states are entirely equivalent in all calculations, 
the large lifetime broadening T causes a destructive interference between the different 
intermediate states leading to the same final state. To understand this, the RIXS 
cross section must be more closely examined. The electron-phonon contribution to 
the intermediate state can be solved by introducing the displaced operators a^ m = 
aim + \/^km/fiwkm> giving a Hamiltonian H = £ km al m a km hu km - A km . Applying 
an ultra-short core- hole lifetime approximation [23, |2H ES] to the intermediate-state 
propagator by expanding the denominator around the resonance energy E res results in 



(f\f\n)(n\f\i) 
z E n -\- E vcs 




with z = hu + Ei — E res + i^. On resonance, where hui + E{ = E res and therefore 
\z\ = T/2, the above expression is only valid when 2\E n — E rcs \/T < 1. At first, this 
does not appear to be justified since the eigenstates of H go to infinity. However, only 
the eigenstates with a finite spectral weight, which are spread in an energy region of 
approximately 2A F around the resonance frequency, are pertinent to the RIXS cross 
section. Limiting the calculations to the lowest order, an effective scattering amplitude 
can be obtained as 



E 



(f\T\n)(n\T\i) 



hu + Ei - E n + 



5(hoj' + E f - hu - E^, (2) 



\(f\f(H - E Tes )f\i) = £ l^/A km ^ m (/|(4 m + a km )|z), (4) 

km 

by omitting the terms that do not give rise to inelastic x-ray scattering. The transition 
operators causing the electronic excitation in the intermediate state effectively cancel 
due to the unilateral consideration of the phonon excitations in the final state. Therefore 
in the limit of a very short core- hole lifetime (the large F limit), only single phonons 
will be excited, and the probability will be proportional to the coupling constant. The 
numerical calculation of the Hamiltonian supports the idea behind the expansion as 
shown in figure |3j For F/huQ = 3 and 5, the lifetime broadening is comparable to the 
width of the x-ray absorption spectral features (inset in figure |3]); the RIXS spectra 
clearly show multiple phonon excitations (figure [3]) and the restriction to lowest order 
in the ultrafast core-hole lifetime expansion is invalid. For F/huo = 10 and 20, the 
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multiphonon features are strongly reduced by the destructive interference, and the 
spectra are dominated by the single-phonon excitations, described by the amplitude 
in Equation (J4]). An alternative way to understand the ultra-short core- hole lifetime 
approximation is by noting that the decay of the core-hole is so fast that the lattice 
has insufficient time to fully respond to the change in charge distribution caused by 
the resonant process. At the K-edge in transition metal compounds, where T > 2 
eV and A/T <C 1, the ultra-short core- hole lifetime approximation is well justified. 
Note that the ratio between the Stokes and anti-Stokes single-phonon features remains 
almost constant for all T. In figure |2j the experimental spectra are compared with 
a numerical calculation where the intermediate states are exactly diagonalized. The 
intermediate-state lifetime broadening T is taken to be 2.0 eV. Some weight has been 
added in the zero-loss peak to account for elastic diffuse scattering. The final spectra are 
further broadened with a Gaussian with a width of 38 meV to account for experimental 
uncertainties, including the instrumental resolution. The relative RIXS intensities of the 
phonon peaks are proportional to Ak m /jWk m (Equation (1)) and therefore are indicative 
of the coupling strength between the phonons and the transient change in charge density 
with respect to the ground state. 

From the intensities of the inelastic phonon peaks, see figure EJ^b), we observe a 
clear resonance behaviour. The largest intensity is observed in the pre-edge region 
(8981 eV), where direct excitations from Is to 3d are made. This causes a change in the 
local 3d charge distribution that directly affects the surrounding ligands. The intensity 
then drops with the main edge (> 8987 eV), which is dominated by excitations into 
the delocalized 4p band that are less effective in coupling to the phonon modes. The 
multiphonon excitations are strongly suppressed due to the short core-hole lifetime. 
The best agreement between theory and experiment is obtained when Ak m for the 
24 and 70 meV features are reduced a factor 0.7 compared to the 41 meV phonon 
excitations. Noting the relative intensities are directly related to the electron-phonon 
coupling strength, the spectra reveals direct information of the relative corresponding 
coupling strengths. High-quality data measured on resonance can therefore provide 
insight into the electron-phonon coupling. The relative ratio between the Stokes and 
anti-Stokes is given by the Boltzmann factor at room temperature, a further confirmation 
that the features are phonon related. 

In conclusion, we have for the first time demonstrated the resonant enhancement of 
phonon excitations at the i^-edge of copper. The experimental features can be explained 
by the coupling between the phonons and the transient change in the charge distribution 
on the site where the resonant scattering process occurs. Due to the fast decay of the 
core-hole, the system has no time to fully respond to the change in charge density, 
and multiple phonon excitations are strongly suppressed. Using an ultra-fast core-hole 
expansion, we demonstrate that the resonant inelastic scattering intensity is directly 
proportional to the electron-phonon coupling strength. Future experiments should 
include high-T c superconductors where knowledge of the relative coupling strengths of 
different phonon modes might shed additional light on their role in superconductivity. 
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We would like to stress that this technique studies the coupling of a specific element 
(Cu in this case) to the lattice. One can also perform experiments at different edges to 
obtain element-specific couplings to the phonon modes. Another aspect not explored in 
this paper is the momentum dependence. The study of the dependence of the intensity 
of the inelastic phonon peaks as a function of transferred momentum would provide 
unique insights into the g-dependence of the electron-phonon coupling. 
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Energy (meV) 

Figure 1. (color online) Resolution curve of the spectrometer. The fit is obtained 
with a linear combination of a Gaussian (dashed) and a Lorentzian (solid-thin), i.e. 
pseudo-voigt function (solid-thick). Full width at half maximum (FWHM) is measured 
to be 38 meV and this number remained unchanged for all incident energies. 
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Figure 2. (color online) (a) Resonant inelastic x-ray scattering spectra at the Cu 
K-edge in CuO for four different incoming energies. The solid lines indicate numerical 
calculations of the RIXS cross section broadened with a Gaussian with a width of 
38 meV. For hw = 8981 eV, the spectra are shown without additional Gaussian 
broadening, (b) x-ray absorption spectrum in fluorescence yield mode at the Cu K- 
edge in CuO. The lines indicate where the RIXS spectra were taken. The length of 
the lines indicates the intensity of the inelastic signal. 
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Figure 3. (color online) Exact calculation of the resonant inelastic x-ray scattering 
for a single phonon mode with energy Hloq = 40 meV coupled to the transient change 
in charge distribution on the site where the absorption takes place (coupling strength 
A is taken to be 0.15 eV). The calculations clearly show the decrease in intensity of 
the multiphonon excitations with decreasing core-hole lifetime (increasing the lifetime 
broadening T). The inset shows the x-ray absorption spectra (XAS). Spectra are given 
for different intermediate-state lifetime broadenings of T/Tlujo —3, 5, 10, and 20. The 
RIXS spectra are normalized to the 40-meV loss feature. The grey lines in the inset 
show a XAS spectrum for T = Twjq/A, clearly showing the phonon excitations. 



